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ABSTRACT 

The evolution with time of the abundance of galaxy clusters is very sensitive to the statistical 
properties of the primordial density perturbations. It can thus be used to probe small devia- 
tions from Gaussianity in the initial conditions. The characterization of such deviations would 
help distinguish between different inflationary scenarios, and provide us with information on 
physical processes which took place in the early Universe. We have found that when the infor- 
mation contained in the galaxy cluster counts is used to reconstruct the dark energy equation 
of state as a function of redshift, assuming erroneously that no primordial non-Gaussianities 
exist, an apparent evolution with time in the effective dark energy equation of state arises, 
characterized by the appearance of a clear discontinuity. 
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1 INTRODUCTION 

One of the most fundamental predictions of the simplest standard, 
single field, slow-roll inflationary cosmology, is that the primordial 
density fluctuations, that seeded the formation of the large-scale 
structure we see today, were nearly Gaussian distributed (see e.g. 
'Creminelli 2003 ; Maldacena 2003 . Lyth & Rodriguez 2005 , S eery] 
|& Lidsey^2005; ,Sefusatti & Komatsu.2007, Such prediction seems 
to be in good agreement with current observations of the cosmic 
microwave background anisotropics (e.g. Slosar et al.||2008) and 



large-scale structure (e.g. Komatsu et al. 



|2011b. Nevertheless, a 



significant, potentially observable level of non-Gaussianity may be 
produced in some inflationary models where any of the conditions 
that give rise to the standard single-field, slow-roll inflation fail. 

The detection of primordial non-Gaussianities would decrease 
considerably the number of viable inflationary models, and it would 
give us an insight on key physical processes that took place in the 
early Universe. Such detection could be achieved through the statis- 
tical characterization of the properties of the large-scale structure, 
namely the bispectrum and/or trispectrum of the galaxy distribution 
(e.g. Sefusatti & Komatsu 2007, Matarrese & Verde 2008. 1, or the 
determination of the evolution with time of the abundance of mas- 
sive collapsed objects such as galaxy clusters (see e.g. |Matarrese et| 
|al.|2000[|Robinson & Baker|2000) . These form at high peaks of the 
density field 6 (x) = Spip and their number density as a function 
of redshift depends on the growth of structure, thus being sensi- 
tive to the dynamics and energy content of the Universe and to the 
statistical properties of the primordial density fluctuations. 

In this work, we address the issue of how primordial non- 
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Gaussianities may affect the determination of the effective dark en- 
ergy equation of state w using the evolution with time of the galaxy 
cluster abundance. Throughout, and unless stated otherwise, we 
consider our fiducial cosmological model to be a flat ACDM model 
with WMAP 7-year cosmological parameters (WMAP-1-BAO-1-//0) 
( [Komatsu et al.|[2011| , namely, a Hubble constant. Ho, equal to 
100/jkms ' Mpc" with h = 0.704, fractional densities of matter 
and baryons today of SI,,, = 0.272, Q.i,h~ = 0.023 respectively, 
a scalar spectral index, equal to 0.963, and we normalize the 
power spectrum so that erg = 0.809. 

2 PRIMORDIAL NON-GAUSSIANITY 

Primordial non-Gaussianity is commonly parametrized by the non- 
linear parameter /(vl and may be written as follows \Lo Verde et al.| 
|2008j), 



Bc(kuk2,h) = {l7lf 



(1) 



where K = k[ +A;2+A:3, ^ is the primordial curvature perturbation and 
^ is an auxiliar function that contains the shape of the bispectrum, 
B(, of Further, cc fc"'"' is the dimensionless power-spectrum, 
while kj = |k,| with k, being the wave vectors. In Fourier space, 
the functions and B( are defined by means of the two and three- 
point correlation functions, 

(2Kf'P,-(k,) 



(^■k,fk,) = (2^)'^«(ki2)- 



47r/tf 



(2) 



(3) 



(fk,fk/k,; = (27TySo(kn3)B((k„k2,k3), 

where k,^ „ = k, + k^ + ... -1- k„. Note that f^i may or may not be 
scale-dependent. In this work only the later case is considered. 
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The bispectrum of ( is the lowest order statistics sensitive 
to non-Gaussian features. Depending on the underlying physical 
mechanism responsible for generating non-Gaussianities, different 
triangular configurations (shapes) will arise. There are broadly four 
classes of triangular shapes or, equivalently, four different bispec- 
trum parametrizations can be found in the literature: Local, Equilat- 
eral, Folded and Orthogonal. Here we shall only consider the first 
two. They are defined as follows: 

- Local shape: It arises from multi-field, inhomogeneous re- 
heating, curvaton and ekpyrotic models. Mathematically, the Local 
shape can be characterized by a simple Taylor expansion around 
the Gaussian curvature perturbation, (q, jSalopek & Bond||1990| 
[Komatsu & Spergel|2001^ , 

3 



C (X) = (a (X) + -/ir' (4 (X) - {& (x)>) 



The WMAP 7-year estimate for /^°™' is ( [Komatsu et al.|2011 1 



rlo 

In. 



32 ±21 (68% C.L.) . 



(4) 



(5) 



The bispectrum for this shape can be derived using Eq. (j4| and it is 
given by ^Lo Verde et al. 2008 , Ko matsu 20 lOj 



(6) 



This quantity is maximized for the so-called squeezed triangle con- 
figuration, i.e. k^ ^ k2 ~ ki. 

- Equilateral shape: It is characteristic of inflationary models 
where scalar fields have a non-canonical kinetic term (for example 
Dirac-Bom-Infield inflation, see |Alishahiha et al.|2004| l. The math- 
ematical expression for the Equilateral shape is i Lo Verde et al.| 
[2008||Komatsu|2010t 

_ 9_ 

"""" " To" 

xl+2(«,,-l)/3 j_ ,^2+(n,,-l)/3j,l+2(,i,,-l)/3;,n,,-l 



^™,„ = —K-^'"'-'^ [-kl (k^hY'-' +perm. 



-2 (kikikY^''"'-"" + fc^+>«..-u/^^.+A«,,-ij/.^»,,-i + pg^^ j ^ (-7) 

This quantity reaches a maximum at k] x k2 ~ k^. Constraints 
from WMAP 7-year set the level of non-Gaussianity for this shape 
at (Komatsu et al.|20Il| l 



/;;f = 26 ± 140 (68% C.L.) 



(8) 



As mentioned before, the abundance of rare objects such as 
galaxy clusters holds relevant information that can be used to probe 
the initial conditions. For such information to be of use, the statis- 
tics of the density perturbation, Sr, smoothed on a scale R, have 
to be characterized. However, we have previously defined the non- 
Gaussianity in the primordial curvature perturbation, rather than 
in the smoothed linear density field. The relation between ( and the 
linear perturbation to the matter density today smoothed on a scale 
R is given by. 



Sr (k, z) = D iz) W(k, R)M{k)T(kX (k, z) , 
where 



(9) 



(10) 



with D(z) being the linear growth factor ( [Komatsu et al. 
2009b, T {k) is the transfer function adopted from Bardeen 
et aLl ( |1986[ l, and W{k,R) is the smoothing top-hat win- 



dow. We use the shape parameter given by |Sugiyama| ( [1995| l, 

F = n,„/iexp [-n,, (i + V2/^/n,„)]. 



Using Eq. ([9]l and the definitions given in Eqs. (|2| and (|3]l, one 
may compute the variance 



J (2n f J (2n f 



r 



dk 



-THk)'P(ik) , 



(11) 



and the three-point function for the smoothed density field ( |Lo| 
[Verde et al.|2008t 



^3) 



(12) 



wiui y = L, \is.i), J i = yy \ 

ki2 = Jk] +kl + 2ki • k2 
Eqs. ([U) and ^ will be 



(^1^2^12)^ 

witii f,- = ^{V.i),Ti= W{ki, R)M(ki)T {k,), K = k, + k2 + k^ and 
k.2 (here a dot represents the scalar product). 
I be of special relevance in the next section, 
since in order to incorporate non-Gaussian initial conditions in the 
prediction of rare objects, one has to derive a non-Gaussian prob- 
ability density function (PDF) for the smoothed density field, 6r. 
This can be done by using a mathematical procedure that enables 
us to construct the PDF from its comulants (see [Lo Verde et aT] 
[20081 [Matarrese & Verde|2Ci08| for details). 



3 THE DARK ENERGY EQUATION OF STATE FROM 
THE GALAXY CLUSTER ABUNDANCE 

3.1 Halo Mass Function 

The comoving number density of virialized halos per unit of vol- 
ume at a given redshift z, dn/dM {z,M), with a mass, M, in the 
range [M, M + dM] is called the mass function. In the presence 
of non-Gaussian initial conditions, the mass function has been es- 
timated using extensions of the Press-Schechter (PS) formalism 
( [Press & Schechter[[1974^ . This formalism asserts that the frac- 
tion of matter ending up in objects of mass M is proportional to 
the probability that the density fluctuations smoothed on the scale 
R = (3M/4np)^ , and above a certain threshold value, d^, can be 
written as. 



dn P d 

{z,M) = -2— 

dM MdM 



dvP (y, M) 



(13) 



where p, <Tm, 5^ and P{v,M) are respectively, the comoving mass 
density, the root square of the variance of the density perturbation 
in spheres of radius R, the critical overdensity in the spherical col- 
lapse model and the PDF of the smoothed density field. The redshift 
dependence of the threshold for spherical collapse as been incorpo- 
rated in 6c (z) = 1.686D (0) (z). For Gaussian initial conditions 
the mass function acquires the following form. 



dn 
dM 



(M,z) ■■ 



2 p 6 AZ) din a-M ^_s2i-y(2al) 



(14) 



n M- (Tm din M 

The cosmological parameters enter Eq. ^14\ essentially through the 
variance and the linear growth factor, as well as, through the critical 
density contrast dc (z). 

There are several prescriptions for chang ing Eq. \14\ in the 
presence of non-Gaussian initial conditions ([Chiu etsi. 1998 



Robinson et al.'2000', |Avelino & Viana|2000| [Fbsalba et al.|2000 
Matarrese et al. 2000}. Here we will adopt that which has been 
proposed by Lo Verde et al.[ ( [2008^ , and which has been shown to 
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(a) (b) (c) 



Figure 1. The number of galaxy clusters per unit of redshift per square degree with mass M > M/„„ = 5 X 10 hT Mq, considering »■ = -1 and different 
levels of non-Gaussianity for (a) Local and (b) Equilateral parametrizations. Panel (c) shows the effect that a change on a constant dark energy equation of 
state parameter, w, has on the number of clusters per unit of redshift per square degree with Gaussian initial conditions {Jnl = 0). 



provide a good fit to results from N-body simulations (see |Wagner| 
|et al.|2010| and references therein), 



dn 
dM 



dlncTu ( 6c iz) SjmCTm 



d\nM 



1 dSiM 

H (Tu 

6d\nM 



cr„ 6 
Sliz) 



(15) 



where 53^ = oc /^^ is the skewness of the smoothed 

density field. If //v/. = 0, then S^m = ^ and Eq. \\5\ reduces to the 
Gaussian mass function. 

Numerical simulations have shown that the PS form of the 
mass function under-predicts the abundance of high-mass objects 
and over- predicts low-mass ones. Therefore, to be in better agree- 
ment with the results of numerical simulations, we follow [Lo Verde] 
|et al.|p008^ and model the departures from Gaussianity using the 
mass function suggested by |Sheth et al.| ( [200T| l. The non-Gaussian 
mass function then becomes. 



drih, 



dM 



{z,MJnl)-- 



dns T dnps /dM(z, M, /a-l) 
dM dnps/dM(z,M,fNL = 0) 



(16) 



The Seth-Tormen mass function has been calibrated using nu- 
merical simulations with Gaussian initial conditions and it is given 
by (,Shethetal.|2001^ , 



diisj 
dM 



(z,M) 



p SAz) 



din M 



(17) 



with a = 0.707, A = 0.322184 and p = 0.3. 

Given Eqs. ^16\ and jl7[ l, one may now compute the num- 
ber of clusters of galaxies per unit of redshift above a certain mass 
threshold, M/;„,, 



dN 



(z,M > Mii,„) 



■ hh-r- (z) 

«Z JMi,„ 



J Ml, 



dM ^iz,MU), 



(18) 



where fsky is the fraction of sky being observed in the survey and 
dV/dz(z) is the comoving volume element which, for a flat cosmol- 
ogy, is given by 



dV 

~d~z 



(z) 



4;r;r(z)?(z), 
dz 



(19) 



with;^' (z) being the comoving radial distance. 



Figures[Ta]and[Tb]illustrate the impact of the primordial non- 
Gaussianities on the number density of galaxy clusters with mass 
threshold M > Mr,,,, = 5 x IO'^/j^'Mq and w = -1, as a func- 
tion of redshift, for the Local and Equilateral parametrizations, re- 
spectively. Figure [Tc| shows how the cluster number density is af- 
fected when a constant dark energy equation of state parameter w 
is changed for the same mass threshold, always assuming Gaus- 
sian initial conditions. The effect of /^.^ and w on the abundance of 
galaxy clusters is quite different. On one hand, increasing w above 
w = - 1 flattens the slope of the cluster abundance above z ~ 0.5, 
which translates in an increase in the number of high-z clusters. On 
the other hand, changes in f^i modify the cluster abundance more 
uniformly in redshift. The difference in behaviour occurs because 
w affects both the volume factor and the mass function, while /wl 
changes only the tail of the distribution of the density perturbations, 
thus modifying just the mass function. 

3.2 Estimation of Wp// 

Figures [Ta| [lb] and [Tc] suggest that the redshift evolution of the 
number density of galaxy clusters in non-Gaussian models could 
be wrongly taken to be the result of an effective dark energy equa- 
tion of state different from the real one, under the assumption of 
Gaussian initial conditions. In order to test this hypothesis, we 
have generated mock catalogues with the expected redshift evolu- 
tion of the cluster number density for different non-Gaussian initial 
conditions in bins of redshift with width Az = 0.1 up to redshift 
2,assuming a mass threshold of M;„„ = 5 x IO'^'/j^'Mq and a nearly 
full sky survey area of 40,000 square degrees. We further consider 
Z^*"' = (-20, 20) and /^f = (-60, 60), respectively for theLocal 
and Equilateral parametrizations mentioned in section [2] The rea- 
son for choosing these specific values for f^t will become clear in 
the next subsection. 

Having generated the mock catalogues with the expected red- 
shift distribution of the number density of galaxy clusters with non- 
Gaussian initial conditions, we then computed an effective dark en- 
ergy equation of state, w^ff, using Eq. l |18[ ) with /jv^ = 0, that mim- 
ics the distribution of the number of galaxy clusters in the presence 
of non-Gaussian initial conditions at the i-th redshift bin, by solving 
the following equation, 



NY 



nzi+Az/2 



dN 
dz 



(z, w = Weff, ffiL =0)dz- 
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Or ' ' ' ' n 1-0 




Figure 2. Figure |2a| shows the value of w^ff that maximizes the abundance of clusters as a function of redshift, with Gaussian initial conditions and M;,„, = 
5 X lO'* /r' Mq (the vertical dotted line corresponds to the redshift z» ~ 0.575 where w» = - 1). Figure [2b| shows the dependence of the of the value of Zt on 
the mass threshold M;;,„. Figures[2c|and[2d|show the reconstructed effective dark energy equation of state w^ff for M > M/,„, = 5 X lO'* /i^'M© and different 
values oi fNL (Local and Equilateral parametrizations, respectively). 



Model 


Local 


Equilateral 


Inl 


-20 20 


-60 60 




0.812 0.806 


0.813 0.807 



Table 1. The computed erg for the Local parametrization and Equilateral parametrization, obtained by demanding the present-day number density of galaxy 
clusters is recovered. 

(Tg = 0.809) is recovered. Table[T|shows the computed erg for dif- 
ferent values of f^t and Local and Equilateral parametrizations. 

Figures [2c] and [2dl show the reconstructed dark energy equa- 
tion of state, Wgff, as a function of the redshift for M/,„, = 5 x 
lO'^/r'Mo and different values of f^t (Local and Equilateral 
parametrizations, respectively). At low redshift, our reconstructed 
Wgff is very close to our fiducial w = -1. But, as we move to- 
wards higher redshifts, our computed w^f f deviates from -1, with 
this effect being more evident for higher values of fm in both 
parametrizations. Further, for f^t > there is a redshift interval 
where w^yy, is undefined, which widens with increasing f^t- For 



-riz,w = -\jNL*Q)dz = Nf° (20) 

Jc,-Az/2 dz 

with respect to w^ff for each redshift bean. Thus, the effective dark 
energy equation of state, Wgff{z), is defined as the value of w which 
reproduces the number of clusters of the mock non-Gaussian cata- 
logue in i-th redshift bin. 

The normalization of the power-spectrum, for models with 
non-Gaussian initial conditions, was done by demanding the 
present-day number density of galaxy clusters in the concordance 
model (ACDM cosmology with Gaussian initial conditions and 
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Mii„ = 5 X IO'^'/z^'Mq this happens in a redshift range centered 
at the redshift z, ~ 0.575, at which the value of w^ff, which we 
will call vv,, that maximizes the cluster abundance is equal to -1 
(see Figure [2a]l. In this interval, there is no Weff that solves Eq. 
\2Q) , since the product of the comoving volume with the integral 
of the mass function, when assuming Gaussian initial conditions, is 
always smaller than the same quantity for non-Gaussian initial con- 
ditions, i.e. fm it 0. If our fiducial cosmological model had a differ- 
ent value for w, then the discontinuity would appear at the redshift 
at which w = w,. On the other hand, for models with /(vl < 0, we 
have also a discontinous w^ff but there are two values of Weff capa- 
ble of reproducing the number of clusters in non-Gaussian models 
in a small redshift interval centered at Zt,- The dependence of the 
value of z, on the mass threshold M;„„ is shown in figure[2b| 

3.3 Estimation of w,// with statistical uncertainties 

The observational estimation of the number density of galaxy clus- 
ters is affected by two sources of statistical uncertainty: the shot- 
noise and the cosmic variance (e.g. Valageas et al. 201 1). The statis- 
tical uncertainty associated with the former increases, for example, 
with the cluster mass threshold, as clusters then become more rare, 
while the statistical uncertainty associated with the later increases, 
for example, as the cosmic volume surveyed gets smaller. However, 
assuming primordial Gaussian density perturbations, and for a clus- 
ter mass threshold above 5 x 10'"*/?"' Mq, it has been shown that the 
magnitude of the contribution of the cosmic variance to the statisti- 
cal uncertainty, in the observed number density of galaxy clusters, 
is always at least an order of magnitude smaller than the contribu- 
tion due to shot-noise, almost independently of the surveyed sky 
area (e.g. see figures 5 and 9 of jValageas et al.|201I^ . 

The existence of statistical uncertainties in the observed num- 
ber density of galaxy clusters will to some extent mask the apparent 
discontinuity on the evolution of w^yy described in the previous sec- 
tion, which could be used to identify the presence of non-Gaussian 
density perturbations. In order to determine the impact of such un- 
certainties, we will re-estimate w^ff in similar fashion to what was 
done in the previous section, but with the inclusion of shot-noise, 
Cn = at the Icr level in each redshift bin. We are able to ne- 
glect the contribution from cosmic variance by setting the cluster 
mass threshold to 5 x lO'^^/i^'M©, and noting that the assumed level 
of non-Gaussianity is relatively small, not affecting much the aver- 
age number density of galaxy clusters with respect to the Gaussian 
case (see Fig. 1). 

Figures |3] and |4] show the Itr statistical uncertainty in the re- 
constructed Wfff, as a function of redshift, when shot-noise is taken 
into account, and as before for a nearly full sky survey area of 
40,000 square degrees. As can be seen, even when including the 
effect of such uncertainty, the apparent discontinuity on the evolu- 
tion of Weff found in the previous subsection is still present (at least 
at the 2cr confidence level) for values of fyt as small as ±20 for the 
Local parametrization, and ±60 for the Equilateral parametrization. 
Clearly, decreasing the survey area would mean that only values of 
/nl with larger modulus could be detected. 

4 CONCLUSION 

In this work, we investigated whether the presence of primordial 
non-Gaussianities has an impact on the estimation of the effec- 
tive dark energy equation of state, when one uses the abundance of 
galaxy clusters as a tool to probe different cosmological scenarios. 
We computed the effective dark energy equation of state, w^ff per 
redshift bin, assuming Gaussian initial conditions, that is capable 



of reproducing the galaxy cluster counts expected in several non- 
Gaussian models, thus constructing a correspondence /jvl w^/f 
for each redshift bin. The most important result of this work is the 
discovery of a redshift interval where no value for the effective 
dark energy equation of state is capable of reproducing the non- 
Gaussian cluster abundance for models with fyi > 0, which is the 
result of there being a w^ff that maximizes the cluster number den- 
sity at each redshift, while for models with f^t < a discontinu- 
ous Wf ff is obtained. The appearance of such type of features may 
thus constitute a new diagnostic of the presence of primordial non- 
Gaussianities. Although, even under the ideal situation where only 
statistical uncertainties are present, their detection is only possible 
for values of /vz. not too far away from those permitted by analy- 
sis of the WMAP 7-year data lI Komatsu et al.pOIl) , it should be 
remembered that /nl may not be invariant with scale. In particular, 
there isn't any particular reason why /;vl could not increase as the 
scale diminishes (see e.g. |Riotto & Sloth|2011^ . 
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Local, fNL = -20, AO = 40000 deg2 




Local, fNL = 20, An = 40000 deg2 




'0.1 0.5 1.0 1.5 1.9 0.1 0.5 1.0 1.5 1.9 

z z 

Figure 3. The reconstructed effective dark energy equation of state for the Local parametiization with statistical observational uncertainties taken into account. 
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Figure 4. The same as in Fig.|3] but for the Equilateral parametrization. 
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